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SPINNING MODE SOUND PROPAGATION IN DUCTS WITH ACOUSTIC TREATMENT 

by Edward J. Rice 


Lewis Research Center 

SUMMARY 

Recent acoustic data have shown larger noise attenuations than predicted for acous- 
tically treated aircraft engine inlets without splitter rings. These encouraging data have 
stimulated a more detailed theoretical study of the acoustic propagation of spinning 
modes in acoustically treated open circular ducts. In addition, the suppressor with 
splitter rings was modeled by using the rectangular approximation to the annular duct. 
The theoretical models were used to determine optimum impedance and maximum atten- 
uation for several spinning lobe numbers from 0 to 50. Some interesting results of the 
analysis are that for circular ducts the maximum possible attenuation and the optimum 
wall impedance are strong functions of the lobe number. For annular ducts the attenua- 
tion and optimum wall impedance are insensitive to the spinning lobe number for well 
cut-on modes. These results help explain why suppressors with splitter rings have been 
quite effective in spite of the lack of detailed information on the noise source modal 
structure. Conversely, effective use of outer wall treatment alone will require expanded 
knowledge of the noise source structure. Approximate solutions are presented to help 
interpret the more exact theoretical results. 


INTRODUCTION 

Noise suppressor design procedures have evolved over the past several years pro- 
viding quite predictable suppressor configurations for aircraft engine inlet and exhaust 
ducts with splitter rings. These procedures have evolved from empirical correlations 
of extensive noise attenuation tests in flow ducts and from theoretical studies of sound 
propagation in ducts with acoustic treatment. The sound propagation theory has used 
quite simple noise input functions such as a plane wave (refs. 1 to 3) or the axisymmet- 
ric least-attenuated mode (ref. 4). These simple models have provided very good design 
procedures as long as splitter rings were used. Not much attention was paid to the sus- 
pected, but largely undefined, complex modal structure of the fan noise source, and it 



really did not seem to make much difference that this complexity was ignored. It has 
been known for several years that spinning duct modes would be produced by rotor - 
stator interaction (ref. 5) and by inlet distortion- rotor interaction (refs. 6 and 7). Yet 
flow duct tests which did not contain this type of mode produced empirical design proce- 
dures that resulted in predictable splitter ring suppressors. 

Recent noise data from aircraft engine inlets showing surprisingly high suppression 
with only outer wall acoustic treatment as well as the strong incentive to reduce noise 
suppressor pressure losses have stimulated renewed interest in the role of modal struc- 
ture in noise suppressor theory. Neither axisymmetric sound propagation theory nor 
flow duct empericism has explained the large noise attenuations which were observed. 

It was hypothesized that the spinning mode character of sound does indeed make a differ- 
ence in the performance of inlet suppressors with treatment on the outer wall alone. 

The propagation of spinning modes in soft-walled ducts has been studied in the past 
(e.g. , refs. 8 to 10 and 11), but a general design procedure has not resulted from these 
studies. It appears that a thorough theoretical and experimental study of spinning mode 
sound propagation in acoustically treated circular ducts is needed. This report repre- 
sents an attempt to initiate this theoretical study to provide rules for optimization of 
duct treatment and to explain some of the observed differences between splitter ring de- 
signs and outer-wall treatment alone configurations. 


SYMBOLS 


A constant multiplier (see eq. (9)) 
a one-half of lined duct height, m 

optimum resistance coefficient for circular ducts (see eq. (36)) 
B optimum reactance coefficient for circular ducts (see eq. (37)) 

A 

c speed of sound, m/ sec 

D circular duct diameter, m 

D m attenuation coefficient (see eq. (60)) 

E acoustic power, W 

f frequency, Hz 

H height between lined walls, m 

J m Bessel function of first kind of order m 
k integer variable 
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L acoustically treated duct length, m 
M uniform steady flow Mach number 
M t lobe pattern tip Mach number 

m lobe pattern number 

P part of acoustic pressure which is function of transverse coordinates (r, $ or 
y, z) (see eq. (2)), N/m 2 

part of acoustic pressure which is function of coordinate between acoustically 

2 

treated surfaces (r or y), N/m 

2 

p acoustic pressure, N/m 

Q 1 + iM(a + iT) 

R amplitude of eigenvalue a 

r radial coordinate, m 

r i inner radius of an annulus, m 

r Q outer radius of an annulus, m 

t time, sec 

\ axial group velocity, m/sec 

6 X 

V axial phase velocity, m/sec 

px 

v acoustic particle velocity (vector), m/sec 

x axial coordinate, m 

Y_ Bessel function of second kind, of order m 
m 

y transverse coordinate (see fig. 3), m 
Z liner acoustic impedance (see eq. (11)), kg/m /sec 
z transverse coordinate (see fig. 3), m 

a complex eigenvalue (a = Re 1 ^) 

/3 transverse dimension (either D or H) 

6 radius ratio, rj/r Q 

£ dimensionless acoustic impedance, Z/pc 

T) frequency parameter, fH/c or fD/c 

9 specific acoustic resistance, real (£) 

O 

p density, kg/m 
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a attenuation coefficient 

t propagation coefficient 

$ angular coordinate, rad 

<p phase of eigenvalue 

X specific acoustic reactance, imag (£) 

u) circular frequency, rad/sec 

Subscripts: 
j j radial mode 

r radial component 

w value at the wall 


INLET SUPPRESSOR DATA COMPARED WITH THEORY 

The motivation for this study came from looking at experimental data on inlet sup- 
pressors without splitter rings. Samples of these data are shown in figures 1 and 2. 

Sound power attenuation spectra from suppressors tested with NASA Lewis fan C 
are shown in figure 1. Fan C is a high tip speed (472 m/sec, 1550 ft/ sec), 1. 74-meter 



Figure 1. - Effect of inlet suppressor configuration on sound 
power attenuation. Fan C, 90 percent speed. 


Rectangular flow duct 
Engine or fan inlet - 
with splitters 
Engine or fan inlet - 
wall treatment only 

Open symbols denote 60 percent fan speed 
Solid symbols denote 90 percent speed 


PO * 



Least attenuated mode theory 

Circular f No flow, 
Rectangular \ m ■ 0 
(annular approximation) 


.4.6.81 2 4 68 10 

Frequency parameter, fp/C 


Figure 2. -Comparison of suppressor theory and data. For rectan- 
gular or annular duct, use p = H; for circular duct, use p = D. 
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( 68 . 3 -in. ) diameter, single fan stage with 26 rotor blades and 60 stator blades (ref. 12 ). 
Two attenuation spectra are shown: wall treatment plus three splitter rings, and wall 
treatment only. Two interesting points about the data should be noted. First, the wall 
treatment alone configuration is very effective. At high frequencies, near the blade 
passage frequency (~2000 Hz), the large attenuation (« 8 dB) obtained cannot be explained 
by axisymmetric sound propagation theory. Second, at low frequencies where multiple 
pure tones dominate («800 Hz peak), not only is the large attenuation unexplainable by 
axisymmetric theory, but the splitter rings make very little difference in the attenuation. 
The insignificant effect of the rings for this type of sound is probably due to the radial 
pressure distribution of the sound. These questions are considered again after the 
spinning mode propagation theory is developed herein. 

Some normalized peak attenuation data points are shown in figure 2 along with 
three theoretical curves. The maximum sound power attenuation is seen to be normal- 
ized by the quantiy L// 3 . The length L is the length of acoustic treatment while /3 
is a radial dimension. The dimension /3 is the duct diameter for open circular ducts 
and the duct height between treated surfaces for annular or rectangular ducts. Note that 
S occurs also in the abscissa (frequency parameter) and should be similarly interpreted 
there. Maximum sound power attenuation will be explained in more detail later (or see 
also ref. 1). For now, it can be simply considered as the normalized peak experimental 


or theoretical sound power attenuation. 


The theoretical maximum attenuation is ob- 


tained with the walls having a particular acoustic impedance (optimum impedance) at a 
particular sound frequency. This optimum impedance may differ for various duct ge- 
ometries or for different modal distributions which will be discussed in more detail 
later. 

Two of the theoretical curves represent the largest possible attenuation for circular 
and rectangular ducts when only the least -attenuated mode is present in these ducts. 

The third curve shows the maximum attenuation if a plane wave were impressed on 
the liner at one end. The least-attenuated mode, often thought of as the most conserva- 
tive indication of liner performance, is seen to potentially damp faster than the plane 
wave input. A combination of modes must be excited to match this plane wave input in 
a soft -walled duct (since a plane wave is not a soft wall solution), and this conglomerate 
is less readily attenuated than the so-called least-attenuated mode. 

It should be noted that the curves of figure 2 do not represent an attenuation spec- 
trum which could be obtained in practice. Each point on the curve represents a liner 
with optimized wall impedance. No real liner could track this optimum impedance over 
a frequency range. At most, a real liner could reach these theoretical limits at one 
frequency only. The theoretical curves were generated without grazing flow (M = 0), but 
they will not change with flow in the region where the attenuation is falling off with 77 
(ref. 3). This lack of dependence on M will be shown in a later section. The L/H = 4 
plane wave solution was shown since it is typical of the splitter ring suppressors tested. 
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Small changes would occur for other values of L/H (ref. 3). 

Some interesting observations can be made by comparing the data with their respec- 
tive theoretical limits. First, the rectangular flow duct data agree quite well with the 
rectangular duct theory. The theoretical curve in this frequency range is not influenced 
very much by the type of modal distribution (plane wave, equal energy, equal amplitude, 
etc. ) existing in the duct. This is due to the small number of propagating modes which 
could exist over this frequency range (i. e. , from two to five modes). 

The second set of data (engine or fan inlet with splitters) might also be compared 
to the rectangular duct theory. A rectangular duct should be a fairly good approxima- 
tion to at least the outer annulus (which should carry most of the acoustic power) of a 
splitter ring suppressor which typically has a radius ratio of about 0. 8. These attenua- 
tion values are seen to be quite low (about one -fourth) with respect to the theory. 

There were several reasons why this difference between theory and experiment may 
have occurred. The first was intentional; the two sides of the annular passage were de- 
signed for different frequencies so as to broaden the bandwidth of sound attenuation. 

Thus, at any given frequency the entire treated area could not be working at optimum as 
was the liner for the theoretical curves. Several uncontrollable phenomena were prob- 
ably also degrading the liner performance. The liner impedance models are not exact, 
and it would indeed be difficult to construct a liner which was truly optimum. Also, the 
boundary layer refraction effect was not considered in the liner design. Whatever the 
reason, it is not too surprising to find experimental attenuation levels well below a 
theoretically predicted maximum. It is logical to consider the theory as an ideal upper 
limit and to use theoretical trends for extrapolation from the observed experimental 
data. This procedure has subsequently produced quite predictable liner designs where 
splitter rings are used. 

The third set of data (wall treatment only) is seen to be an order of magnitude above 
the axisymmetric propagation theoretical estimate for circular ducts. It is probable that 
data could be obtained which would be even higher if the walls had been truly optimized. 
These data are quite easily explained if one accepts that the sound is mainly concentrated 
in spinning modes which damp faster than the axisymmetric mode theory of figure 2 
would suggest. However, spinning modes should also exist with splitter rings. At first 
one might expect that these modes should also damp faster than the axisymmetric modes. 
But the data of course deny this conclusion. It will be shown after the development of the 
spinning mode theory that the spinning mode properties are largely irrelevant in the 
splitter ring case, but are the crucial factor in the case of wall treatment only. 


SPINNING MODE PROPAGATION THEORY 
Certain assumptions exist for the theoretical development which follows. The 
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steady flow is uniform everywhere in the duct. Boundary layers are sufficiently small 
so that their effect may be neglected. The soft walls are assumed point reacting and of 
uniform impedance. These assumptions imply that no waves are propagating within the 
liner backing cavity (prevented by honeycomb cells) and that any perforations are closely 
spaced compared to a wavelength of sound. No end reflections are considered; that is, 
the soft duct is sufficiently long to render the reflections at the duct termination unim- 
portant, or the mode is well enough "cut-on" so that no end reflection will occur. A 
cut-on mode is defined here as a propagating mode. The frequency parameter r\ is 
sufficiently high so that axial propagation occurs. Only very near cutoff will this as- 
sumption (no end reflection) be a problem. Below cutoff the modes carry very little 
energy due to their out-of-phase pressure-velocity relation, and well above cutoff the 
end impedance (ref. 13) is such that very little reflection will occur. It is further as- 
sumed that the sound pressure level is sufficiently low so that the linearized acoustic 
wave equation is valid. 

The acoustic wave equation is expressed as 

v 2 p=J- + 3fE + MJ?IL + M 2 ^ (1) 

c 2 at 2 c dx at ax 2 

p 

TVm nranmofri r* nnfofimi iico^ io ofmwn f i orii r»o *3 Tr» fz-xllrvuTir-irr nnrnn-vnnKc fVio V7^ 

w* V»*W*V MUV' v% »» li jl.±± V V* JU1 Vi IV XVAAV M JLi y ViiV * 




(b) Rectangular. 

Figure 3. - Geometry for analytical models. 
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operator is defined for Cartesian and cylindrical coordinate systems. If separation of 
variables solutions are sought in the form of 


P = 


QQ OO 

Ivl 


3=1 3=1 


iwt-(w/c)(a.+iT.)x 
p . e 1 J 


(2) 


where is a function of the transverse dimensions (y, z or r,4>)then, from the line- 
arized momentum equations, the acoustic particle velocity can be expressed as 


V j 


i Vp. 

pcuQj 


(3) 


where 


Qj = 1 + iM(a. + iTj) 


(4) 


The j subscript is used to denote the infinite number of possible modal solutions. 


Circular Duct 

For a circular duct (cylindrical coordinates), the left side of equation (1) can be ex- 
pressed as 


v 2 p. Ii.fr 

r 3r V 3r/ r 2 ; ^2 Sj[ 2 

The remaining term in equation (2) is then sought in the form 

Pj(r, $) =^’ j (r)e" im4. 

Thus, using equations (1), (2), (5), and (6) results in 



(5) 


( 6 ) 


(7) 
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where 



- M 2 )(a. + ir.) 2 + i2M(o. + ir.) + l] 

J J J J 


( 8 ) 


is the equation that links the eigenvalues with the axial attenuation and propagation 
coefficients (o, r). Equation (7) is recognized as Bessel's differential equation with so- 
lution (ref. 14) 


o^j(r) = J 


m 



+ 


A i Y m 

i m 



(9) 


where the a j and Aj must be determined by radial boundary conditions. A second 
constant of integration which usually appears before J m has been set equal to 1. 

Equation (9) is valid for annular or open circular ducts. However, for added sim- 
plicity in this report, the annular duct case (splitter rings) will be approximated by a 
rectangular duct development in the next section. Of immediate interest, however, is 
the case of an open circular duct. For this geometry, since Y m (^r/r Q ) becomes in- 
finite as r - 0, the condition must be imposed that 


Aj = 0 (all j) (10) 

The remaining radial boundary condition involves the wall impedance at r = r Q . The 
acoustic impedance can be expressed as 

( 11 ) 


where v is the radial velocity within the wall (axial velocity within the wall is as- 
sumed to be zero). When the continuity of radial particle displacement is used the ra- 
dial particle velocity just outside the wall for the j** 1 radial mode is (ref. 15, p. 314, 
eq. (34)) 


Z =-£- 
v rw 


v rj - S'™-) (12) 

Note that particle velocity continuity is not satisfied unless Qj = 1 (M = 0). This is due 
to the artificial slip condition of the fluid at the wall. Using equations (2), (3), (6), (11), 
and (12) and defining a new normalized acoustic impedance give 
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(13) 


C 


_z_ 

pc 





where 



(14) 


Equation (13) becomes 


; - (1 

cv-J i(at-) - mJ (a.) 
j m+rj' m v y 

when equations (9) and (10) are used. From this equation (15) the eigenvalues o-j can 
be determined for a given wall impedance, lobe number m, steady flow Mach number, 
duct size, and frequency. If continuity of particle velocity had been used, would 
have appeared to the first power in equation (15). 

The remaining necessary boundary conditions have been subtly satisifed due to the 
choice of solution forms in equations (2) and (6). Tangential continuity of pressure and 
velocity are satisfied if m is an integer in equation (6). Only positive x traveling 
waves are considered by the choice of the x function in equation (2); end reflections 
are neglected. 


Rectangular Duct 

Due to the relative simplicity of the sinusoidal functions compared to the Bessel 
functions, the annular duct (splitter ring case) will be approximated by a semi-infinite 
rectangular duct. 

For a rectangular coordinate system the left side of equation (1) becomes 


v 2 p = ^P + ^E + 3!e 
3x 2 3y 2 3z 2 


(16) 


The analogy to equation (6) is 
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(17) 


Pj(y> z) =^j(y)e 


imz/r 

1 o 


where z/r Q has taken the place of <f>. The spinning wave pattern repeats continuously 
over z distances equal to the outer circumference of the annulus. 

As previously mentioned, this rectangular duct solution should be a good approxima- 
tion to the true annular solution if the radius of curvature is large and if the radius ratio 
is sufficiently close to 1. The conditions are usually satisfied at least for the outer pas- 
sage of a splitter ring inlet which passage might be expected to carry most of the acous- 
tic power. 

Using equations (16), (17), and (2) in equation (1) results in 


d 2 ^ /a \ 2 

i + (-1 ] ^ = 0 (18) 

dy 2 W J 


where 



(1 - M 2 )(a. + ir-) 2 + i2M(cr. + ir.) + 1 

J J J J 



(19) 


which links the eigenvalue cVj to the axial attenuation and propagation coefficients. 
Equation (18) has solutions 


^j(y) = cos 




( 20 ) 


which represent symmetric and asymmetric modes, respectively. Following a proce- 
dure analogous to equations (11) to (15) provides the wall impedance equations: 


C = 



cos Oj 


sin Oj 


for symmetric modes, and 


S = 


-i^Q 2 sin Oj 


a- cos o- 

J J 


( 21 ) 


( 22 ) 


11 


for asymmetric modes where the impedance condition was applied to y = ±a and where 
now 


ua 
T) = 


nc 


fH 

c 


(23) 


with H being the duct height between the two walls with acoustic treatment. Again, 
either equation (21) or (22) will provide the eigenvalues when £, m, M, H, and f 
are provided. 


Acoustic Power Attenuation 


In the more general development where multimodal propagation is considered, the 
acoustic power attenuation is quite complicated (refs. 1 and 3). The complete pressure 
and the three velocity components must be constructed from the combination of all of 
the modes. The modal participation is determined by the input waveform from the noise 
source. The acoustic intensity is then determined and, since products of pressure and 
velocity are involved, modal cross coupling must be considered. The intensity must be 
integrated over the cross section at two axial stations, say x = 0 and L, since the ra- 
dial distribution of intensity may change with axial distance. Finally these integrated 
intensities which form the acoustic power can be used to determine an attenuation from 


A dB = 10 



(24) 


For the case of a single mode the situation is much simplier. Since a mode does not 
change relative shape as it travels down the duct, no integration is necessary. The pro- 
cedure just outlined reduces to 


A dB = 10 


log 


-2a)Oj L/c 


) 


coa.L 

- 8 . 68 — 1 — 
c 


(25) 


or 


— — — = -17. 4 TtvjO- (26) 

L/D 3 


for circular ducts and 
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A dB yr, A 
= “17. 4 Jiria. 

L/H J 


( 27 ) 


for rectangular ducts. Although equations (26) and (27) look the same, the q terms are 
defined differently, and cr is different for rectangular and circular ducts. 

The damping cr and propagation Tj coefficients can be determined from equations 
(8) and (19) to be 


°J + lT J = 


-iM + i '1/1 - (1 - M ) 


1 - M 


V3 




for circular ducts and 


v l v 



(28) 


(29) 


for rectangular ducts. The last term under the radical in equation (29) was obtained 
from the expression 


M. = — - = —^2 — (30) 

1 me m(l - 6) 

which defines the tip Mach number of the spinning lobe pattern where 6 is the equivalent 
annular duct radius ratio. 

The calculation procedure and the results of the previous theory for a single mode 
are now presented followed by a return to the theory to obtain approximate solutions 
which will help explain the exact solutions. 


CALCULATION PROCEDURE 

The calculation procedure for cylindrical ducts is outlined here. The rectangular 
duct procedure would follow a similar pattern except for the added simplification of 
having only sine and cosine functions. The method used to handle the Bessel functions 
with complex arguments is also outlined. 
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Calculation of Maximum Attenuation and Optimum Impedance 


The main complications center on equation (15) which relates the infinite number of 
eigenvalues Oj to the normalized wall acoustic impedance £. If a multimodal analysis 
were made (ref. 1) the impedance could be specified and the eigenvalues determined by 
iteration, series expansion, etc. For a single mode an inverse procedure greatly sim- 
plifies the problem. The eigenvalue is specified and the impedance can be immediately 
calculated. Since a major point of this report is to determine the optimum impedance, 
this procedure was used. Considerable care was taken to insure that the selected eigen- 
value belonged to the particular mode of interest. 

For the actual calculations the starting point was one step removed from equa- 
tion (15). Since contours of constant attenuation were desired and A dB is proportional 
to 0^ (eq. (26)), a value of cr was first assumed. Values of Tj were then chosen 
which now define the complex eigenvalue <*j through equation (8) or (28) which are ver- 
sions of the same expression. The eigenvalue is then inserted in equation (15) and the 
complex acoustic impedance is determined. The values of M, tj, and m are held fixed 
for the calculation. These could later be varied to perform a complete parametric 
study. 

A sample of the calculation results can be seen in figure 4 where lines of constant 
Oj and Tj are shown on the wall impedance plane. Note that both components of imped- 
ance have been divided by the frequency parameter tj. This makes the plots similar for 
any tj provided M = 0 or that the mode is well cut-on if M * 0. This last point will 



Reactance, xlV 

Figure 4. - Sample results showing optimum impedance 
definition. Rectangular duct: uniform steady flow 
Mach number, M * 0; lobe pattern number, m = 0-, 
frequency parameter, y ■ fH/C = 5. 
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become more apparent in a later discussion. The initial calculation used a very small 
ON and a carefully selected Tj to provide a large resistance and near -zero reactance 
(hard wall eigenvalues well known) to insure that the calculation was proceeding in the 
proper mode. The choice of Tj determines the mode for a fixed value of cr. The cal- 
culation then proceeds by increasing the value of a until the teardrop-shaped constant 
damping curves degenerate to a point. This point then represents the optimum imped- 
ance for this mode, and the value of Oj determines the maximum possible sound power 
attenuation (from eq. (26)) for the given value of 77 . Actually, only one value of 77 need 
be considered for the cases of most interest since the eigenvalue at the optimum will 
not change with 77 . This fixed Oj can then be used to determine the maximum damping 
at any other 77 by using equations (28) and (26). The cases of most interest again are 
the M = 0 and that for M * 0 and the mode is cut-on. This really means that Qj in 
equation (15) and defined by equation (4) is not a function of 77 but can be approximated 
by 


(31) 

1 1 + M 

for well cut-on modes. From equation (4), Qj = 1 for M = 0 at any value of 77 . This 
is derived in the discussion preceding equation (54). 

If is pushed above the maximum defined previously, one finds that the results 
generate a new set of contours collapsing to a point at some lower impedance; these con- 
tours represent the next higher radial mode. 

This procedure provides the maximum damping and optimum impedance for the least 
attenuated mode if the starting values of o, and r- are properly chosen. 

J J 


Calculation of Complex Bessel Functions 

The Bessel functions with complex arguments could be calculated by the procedures 
given in reference 1. However, for the multitude of calculations required in the previous 
section it is most efficient to inspect one calculation before proceeding with the next. 

This requires the rapid communication which can conveniently be obtained from a desk 
calculator. The desk calculator often cannot handle the large subroutines required by 
the methods of reference 1. Thus, the simplified procedures which follow were devel- 
oped so that the Bessel functions themselves need not be calculated; only the ratio of two 
Bessel functions need be considered. 

Equation (15) may be rewritten as 
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(32) 


S = 


invtf 


a 


^aS a j^ 


m 


The ratio J m+ i/J m can be calculated as a single unit without ever considering each 
part. The Bessel function recursion relation can be written (ref. 14) as 


2k 


+ J k+l(®p J k^j^ 

J J a J 


(33) 


or 


_ 2k _ J k+l^ Q, p ( 3 , 

J k (Q, j) J k (Q 'j ) 

For any ffj, a sufficiently large value of k can be chosen so that the series represen- 
tation of the Bessel function is rapidly convergent and the last ratio in equation (34) can 
be expressed as 


J k+l (cy j ) w tt J 

J k (ffj) 2(k + 1) 


1 +• 



(k + l)(k + 2) 


(35) 


If the first term of equation (35) is considered with large k it is seen to be approxi- 
mately the inverse of and quite small compared to the first term of equation (34). In 
fact, the second term of equation (34) could be neglected to start the recursion calcula- 
tion although for this report equation (35) was used. Equation (34) is repeatedly used 
with decreasing integers k until k = m + 1. Note that at each step the left side pro- 
vides the inverse of the last term of equation (34) needed for the next step. Finally, 
equation (34) provides the inverse of the Bessel function ratio needed in equation (32). 
This procedure provides excellent results with k * 4 1 o-j | . 
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RESULTS OF SPINNING MODE PROPAGATION THEORY 


The following sections present some results of the previously discussed sound prop- 
agation theory. The discussion pertains to the least attenuated mode for a given spin- 
ning lobe number. 


Circular Duct 

In figure 5 the maximum possible attenuation for the least attenuated mode is shown 
as a function of the frequency parameter r\. Several curves of constant lobe number 
(m = 0 to 50) are shown. Several loci of constant lobe pattern tip Mach numbers 
(Mj. = 1 to 4) are also shown superimposed. This tip Mach number (M^ = n^/m) can be 
easily derived from lobe number, geometry, and frequency considerations. The inter- 
esting point to note in figure 5 is that the maximum possible attenuation (peak attenuation 
for an optimized system) is a very strong function of the pattern lobe number m with 
larger values of m having larger values of attenuation. This may be the explanation 
for the apparently high experimental sound attenuations shown in figure 2 for circular 
ducts with wall treatment only. These data were high with respect to the axisymmetric 
theory (m = 0) but not necessarily high with respect to the spinning lobe theory. 



Figure 5. - Maximum attenuation for least attenuated 
spinning modes in circular duct. Uniform steady flow 
Mach number, M = 0. 
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It should be cautioned that a comparison of the data from figure 2 with the theory of 
figure 5 probably does not determine the spinning lobe number because the experimental 
suppressors were not necessarily optimized. Much higher experimental attenuations 
might have been attained with some other liner wall constructions. 

The optimum wall impedance which goes along with the attenuations of figure 5 is 
given by 


BqV 

Rc 

(1 + M) 2 

nt] 

A” 


(1 + M) 2 

nrj 


(36) 


(37) 


with the coefficients Bg and given in table I for various values of lobe number and 
where R is obtained from o'=Re 1< ^. 

The corresponding eigenvalues at the optimum impedance are given in table II. Note 
that as the lobe number increases the optimum acoustic resistance (see Bg) decreases 
and the optimum acoustic reactance (see B^) moves toward zero (the liner tuned fre- 
quency). This would imply that with spinning modes the peak attenuation would occur 
at a higher frequency (for a given liner construction) than would the axisymmetric mode. 
The approximations given in equations (36) and (37) are for well-cut-on modes with 

— < 1 (38) 

H7] 

where R is the magnitude of the eigenvalue for the mode. The locus of R/jrrj = 1 is 
shown in figure 5. The correction for steady flow Mach number will become apparent 
after a later discussion of approximate results. 

It is apparent from the results shown in figure 5, by equations (36) and (37), and in 
table I that, for optimum design of circular duct acoustic suppression, much informa- 
tion must be known about the noise source. As a minimum, at least the lobe number of 
the dominant mode must be known. Preferably the distribution of acoustic power within 
the radial modes should also be known, since radial modes other than the least attenu- 
ated mode would cause changes to be made in the numerical results of figure 5 and ta- 
ble I. The general conclusion that sound attenuation depends on lobe number would not 
be changed for radial modes other than the least attenuated mode. 
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Rectangular Approximation to Annular Ducts 


The results of the rectangular duct theory, intended to approximate the annular duct 
theory, are shown in figure 6 . Again the maximum possible attenuation normalized by 
the duct L/H is shown as a function of frequency parameter 77 = fH/c. Several lobe 
tip Mach numbers and lobe numbers are shown. The lobe tip Mach number Mj., fre- 
quency parameter 77 , and lobe number m relation shown in the figure comes from 
equation (30) with radius ratio 6 = r^/r 0 = 0. 8 . 



Figure 6. - Maximum attenuation for least attenuated 
spinning modes for rectangular approximation to an- 
nular duct. Radius ratio, 6 * rj/r 0 - 0. 8. 


In contrast to figure 5 for circular ducts, the rectangular duct theory shows that for 
moderately cut-on modes (M t > 1. 4) the spinning lobe theory collapses rapidly to the 
axisymmetric propagation theory (m = 0). Another way of viewing figure 6 is that for the 
usual range of interest for splitter ring suppressors (1 < 77 < 2 ) little difference in at- 
tenuation is observed between 20 lobes and 0 lobes. Thus, the attenuation obtained from 
splitter ring suppressors is relatively insensitive to spinning lobe number if the modes 
are well cut-on. This conclusion is supported by the recent work of Yurkovich (ref. 11). 

The optimum wall acoustic impedance for the least attenuated spinning mode in a 
rectangular duct can be given by a single set of equations: 
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0.93 


(39) 


(1 + M)' 


_R 

nr) 


< 1 


-0.74 


(1 + M)‘ 


2L *< i 


nr) 


(40) 


which are applicable for any lobe number m. These equations again represent results 
for the least-attenuated mode. The optimum eigenvalue was given by R« 2.364 and 

(p « 28°. 

Thus, if a splitter ring inlet were being designed for the least-attenuated spinning 
mode and if this mode were well cut-on, the optimum liner design would be insensitive 
to the spinning lobe number as would be the attenuation. This interesting conclusion 
helps to explain the experimental results for the splitter-ring inlets shown in figure 2. 
Although the noise may be distributed in spinning lobe patterns, the axisymmetric 
(m = 0) theoretical attenuation may be close to an upper limit for this type of suppressor. 

The insensitivity of splitter -ring suppressor design and attenuation to spinning lobe 
number may also explain the quite good results obtained from previous suppressors in 
spite of the meager information available about the noise source. 


Comments on Rectangular Approximation to Annular Duct 


The use of the rectangular wave equation (eqs. (1) and (16)) to replace the more 
exact cylindrical coordinate system of the annular splitter ring passage should create 
no problem for the limited scope solutions of interest in this report. For the high ra- 
dius ratios of interest, curvature effects should be negligible. In reference 15 radius 
ratios from 0 to about 0. 8 were used in a cylindrical system with constant frequency 
parameter r), and no change was found in the propagation or attenuation coefficients 
when > 0. 4. This implies that curvature effects were negligible above this ra- 

dius ratio 0. 4. In reference 15 only nonspinning mode results were shown, and it is 
conceivable that spinning modes with high spiral angles may be more sensitive to curva- 
ture. However, at least for well cut-on spinning modes, there is very little spiral angle 
in a cylindrical duct which can be seen by looking at the axial phase velocity 


v 


px 


c(l + M) 


r 1 + i + M /R 

\ 2 

■\ cos 2cp 

(41) 

2 \7tt; 

!/ J 



or the axial group velocity 
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(42) 



c(l + M) 



which are slightly greater and less, respectively, than those of a plane wave (no spiral 
angle) in a hard duct which are both equal to c(l + M). This is true since R/nr] is 
small for a well cut-on mode. A recent study (ref. 11) has shown results for spinning 
modes in an annular duct. For high radius ratios the annular duct eigenvalues were 
shown to converge to the rectangular duct eigenvalues for hard wall ducts. 

Another problem does arise, however, with regard to what modal structure of 
sound is excited within the annular liner passages. The results of this report are for 
the least attenuated mode in either the circular or annular duct. This represents the 
most conservative single mode result. The problem in the splitter ring passage can be 
illustrated in figure 7. The first radial mode (cylindrical coordinates) with small and 



Figure 7. - Pressure profiles as function of circumferential lobe 
number m. 


large lobe numbers are sketched as input to the splitter ring passage. For small enough 
lobe number, the pressure is more uniform over the annulus and the least attenuated 
mode will dominate in the carrying of sound power. The results of this report should be 
valid for this case. As lobe number is increased, more higher radial mode content will 
be excited in the annulus, and the least attenuated mode approximation will underesti- 
mate the true attenuation because these higher radial modes have higher damping coef- 
ficients. 

If the lobe number is increased even further (dashed line, fig. 7), the pressure 
amplitude at the splitter ring approaches zero and the ring is completely ignored in the 
sound propagation. For this case the open circular duct (no splitter rings) solution be- 
comes valid. This can be further emphasized by repeating equation (9) 
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(9) 


^j(r) = J 


m 



+ 



which represents an annular duct solution, and by observing the value of (first ra- 
dial mode) in figure 8. Figure 8 results are from reference 16 for hard wall annular 



Figure 8. - Coefficient of Bessel function Y m 
related to lobe number and radius ratio for 
annular ducts with hard walls. 


ducts. For a given radius ratio 6, as the lobe number m is increased, Aj rapidly 
drops toward zero. Of course as A^ — 0, the open circular duct solution is approached, 
that is, 

(43) 

and the presence of the splitter ring becomes irrelevant to the solution. The same phe- 
nomenon occurs for higher radial modes but at higher lobe numbers. This type of modal 
pressure distribution may account for the noise suppressor data at low frequencies 
shown in figure 1. Recall that for frequencies of about 1000 hertz and lower the splitter 
rings produced very little added attenuation. 
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INTERPRETATION OF THEORETICAL RESULTS 


USING APPROXIMATE EXPRESSIONS 

It is useful to derive approximate expressions from the preceding theoretical anal- 
ysis to clarify the basic difference between noise suppressors with splitter rings and 
with wall treatment alone. Approximations stated earlier without proof are derived 
here. The approximations are mainly for well cut-on modes, which are of the most in- 
terest, since they carry the bulk of the acoustic power. 


Well Cut -on Modes 

Some preliminary approximations are required; the first being applied to 


o + ir = 



1 - M* 


(29) 


for the damping and propagation coefficients for rectangular ducts. The subscripts are 
dropped and it is understood that the least attenuated mode is considered. For a well 
cut-on mode, 


(1 - M 2 ) 




(44) 


which allows the approximate expression 

^1 - x « 1 - - x « 1 (45) 

2 


to be used. Also, 


a 2 = R 2 e i2 ^ = R 2 (cos 2cp + i sin 2 cp) 
Therefore, equation (29) can be written as 


(46) 
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(47) 


o + ir * 



sin 2cp + i ■ 


1 + M 


1 

2 



cos 2 <p + 


M t 2 j 


The damping coefficient is thus 


a 


- (— \ sin 2 (p 

2 W 


and the propagation coefficient is 


T » 


1 

1 + M 



(48) 


(49) 


For circular ducts, equation (28) is seen to be the same as equation (29) except that 
the Mf term is missing and that the lobe pattern tip Mach number information is in the 
( a/nrir term. For circular ducts, the equations can be immediately determined from 
equations (48) and (49) by neglecting the term: 


a 



sin 2cp 


(50) 


T ~ 


1 

1 + M 


A (*)* cos 2<p 

2 \n-nJ 


(51) 


Note that the Mach number M has cancelled out of the damping equations (48) and (50). 

It will be shown that the eigenvalue a = Re 1 ^ is constant at the optimum point (litniting 
contour in fig. 4) even with a steady-flow Mach number. Thus, the maximum possible 
attenuation is not a function of Mach number but depends only on the frequency param- 
eter r). 

Another approximation required is for the Qj term in the wall impedance equations. 
This was given in equation (4) as (again dropping subscripts) 


Q = 1 + iM(a + ir) 


(52) 


Note that in equation (47) the eigenvalue had to be carried to second power to get the 
damping coefficient a. However, for well cut-on modes equation (47) is dominated by 


o + ir « 

1 + M 


( 53 ) 
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Thus, equation (52) can be estimated by 


Q* — — -5- < 1 

1 + M 7T77 


(54) 


Equation (54) can now be used in the wall impedance equations (21) and (15) to provide 


C(1 + M ) 2 

V 


iff cos a 
a sin a 


(55) 


for rectangular ducts and 


C(1 + M ) 2 

V 


iffJ m (a) 




(56) 


for circular ducts. 

Rectangular approximation to annular ducts . - Some interesting conclusions can be 
drawn from equation (55) for rectangular ducts. Consider the left side as an effective 
impedance and note that the right side is a function only of the eigenvalue. Thus, an 
optimization procedure as explained in the CAUJ ULATION PROCEDURE section and as 
shown in figure 4 need be done only once if the coordinates in figure 4 are considered as 
effective impedance instead of as impedance only. The eigenvalue a at which the con- 
stant a contours collapse to a point (optimum effective impedance) will always be the 
same. Since a - Re 1 ^ is always constant at the optimum effective impedance, the con- 
clusion stated after equation (51) is now justified. Since the damping coefficient (a in 
eq. (48)) is not a function of Mach number M and a is constant, the optimum damping 
is a function only of the frequency parameter rj. To carry this conclusion one step 
further, from equation (27) the normalized maximum possible attenuation is a function 
only of the frequency parameter 77 for well cut-on modes and is given by 

A dB ^ 12. 8 

L/H 77 

where the optimum rectangular eigenvalue R = 2. 364, cp = 28° was used. 

Going back to equation (55), if the optimization procedure is performed at M = 0 
and 77 = 1 , it is found that 


C = 0. 93 - 0. 74i 


(58) 


Thus, in general, 
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( 59 ) 


^ * (0.93 - 0. 74i)77 
(1 + M) 2 


which justifies equations (39) and (40). 

Recall that the lobe number m did not occur either in the damping coefficient 
(eq. (48)) or in the wall impedance (eq. (55)). Thus, for well cut-on modes the maxi- 
mum possible attenuation and the optimum wall impedance are not functions of the spin- 
ning lobe number, a result which was also shown by the more exact theoretical results. 

Circular ducts. - Some important differences between circular and rectangular 
ducts can be deduced from the approximate equations. For circular ducts the right side 
of the effective impedance equation (56) is seen to be a function of the lobe number m. 
Thus, although the optimization procedure of figure 4 can be performed for M = 0 and 
rj = 1 and then applied to all M and 77 , the procedure must be repeated for each value 
of m of interest. This results in the optimum impedance coefficients of table I and the 
optimum eigenvalues of table II being functions of lobe number. This also implies 
through equations (50) and (26) that the normalized maximum possible attenuation 
A dB(L/D) is a function of lobe number. It is still not a function of steady flow Mach 
number for R/irq < 1 in figure 5. Equations (26) and (50) can be combined to yield 


AdB_ D m 

L/D 77 


(60) 


where D m is an attenuation coefficient dependent on lobe number m and given in ta- 
ble III. 

Note that equation (60) used in conjunction with table III yields an approximation for 
large values of the frequency parameter 77 . This inverse function is valid for the por- 
tions of the curves in figure 5 which are straight lines. 


Cutoff Modes 

When a mode is cut off it is of less interest since it carries no power and is atten- 
uated in even a hard walled duct. For completeness of the approximate equations, the 
attenuation of the cutoff modes will be developed. The (c^Atj) 2 term in equation (28) 
or (29) will dominate the radical and with M = 0 the result will be 

a « R_cosj£ ( 61 ) 

7T77 


for either rectangular or circular ducts. 
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Using equation (26) or (27) yields 


-17. 4R cos <p (62) 

L/ ft 


where /3 = D for circular ducts and /3 = H for rectangular or annular ducts. Note that 
this expression does not contain the frequency parameter rj and represents the hori- 
zontal portions of the curves in figures 5 and 6. This behavior has been observed before 
in references 1, 2, and 17. For M = 0 the same eigenvalues given following equa- 
tion (40) for rectangular ducts and in table II for circular ducts can be used. However 
for M * 0, the situation is much more complicated. The simple expression (eq. (54)) 
for Q in the wall impedance equation is no longer valid: Q will be a function of rj, and 
instead of just a simple stretching or compression of the equal damping curves in the 
impedance plane (see fig. 4) a rotation will also occur. The eigenvalues will thus be 
changed from those of a well cut-on mode. 

An example of the effect of Mach number on the maximum possible attenuation for 
the least attenuated mode in rectangular ducts can be found in figure HI- 13 (p. 87) of 
reference 4. It should be noted that only at low frequency parameter values is this max- 
imum attenuation a function of Mach number where attenuation increases for inlets and 
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previously discussed. 


CONCLUDING REMARKS 

Suppressors with outer wall treatment alone have been observed to provide much 
more attenuation than previously expected. Such expectations were based on the extrap- 
olation from splitter ring suppressor measurements or least attenuated axisymmetric, 
circular duct theory. The theoretical analysis of the least attenuated spinning mode with 
sound absorbing walls has revealed that the experimental results could have been antic- 
ipated qualitatively. Approximate expressions for the attenuation of well cut-on modes 
were presented to illustrate the theoretical results. To further pursue this approach 
more information is required about the spinning modal structure of the sound emanating 
from the source. 

It was consistent with the spirit of this report, then, to make all necessary assump- 
tions required to point out as simply as possible the gross behavior of both splitter ring 
and wall only suppressors. Many of the conclusions made here could possibly be chal- 
lenged by postulating modes near cut-off or more complex radial mode structures. 
However, the experimental data do provide credibility, if not precision, to the assump- 
tions made here. 
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CONCLUSIONS 


Within the restrictions of the assumptions made within this report the following con- 
clusions were made: 

1. Noise suppressors with splitter rings are relatively insensitive to the spinning 
mode lobe number. This explains why quite simple nonspinning mode theory and rec- 
tangular flow duct tests have provided quite predictable splitter ring designs in the past 
even where spinning modes were possible. 

2. Circular noise suppressors with wall treatment only are very sensitive to spin- 
ning mode lobe number. This explains why axisymmetric cylindrical duct theory does 
not agree well with experimental data taken with fan noise sources. 

3. Due to conclusion 2, much more detailed information about the fan noise source 
is required to best take advantage of the great potential of wall-treatment-only sup- 
pressors. 

4. Approximate expressions for sound attenuation for well cut-on modes have been 
presented which provide insight into the difference between the two types of suppressors. 

5. In a splitter ring suppressor (or an annular exhaust duct), the presence of the 
inner wall can under some circumstances be completely ignored by the propagating 
spinning mode. This occurs when the radial pressure distribution is skewed toward the 
outer wall and very little pressure oscillation occurs on the inner wall. A high spinning 
lobe number and low radial mode cause this to occur. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 4, 1975, 

505-03. 
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TABLE I. - OPTIMUM RESISTANCE AND REACTANCE COEFFICIENTS 


FOR SPINNING MODES IN CIRCULAR DUCTS 



Lobe pattern number, m 

0 

1 

3 

5 

7 

10 

20 

50 

Optimum resistance coefficient, B 0 

0.890 

0. 644 

0.448 

0. 357 

0.303 

0.251 

0. 170 

0.097 

Optimum reactance coefficient, 

0. 38 

0. 221 

0. 121 

0. 095 

0.076 

0.060 

0. 038 

0. 020 


TABLE II. - EIGENVALUES AT OPTIMUM IMPEDANCE FOR 


SPINNING MODES IN CIRCULAR DUCTS 




Lobe pattern number, m 



0 

1 

3 

5 

7 

10 

20 

50 

Amplitude of eigenvalue, 

R 

3. 278 

4.691 

7. 299 

9. 759 

12. 088 

15. 547 

26.662 

58.620 

Phase of eigenvalue, <p, 

deg 

23. 19 

18. 07 

13.49 

10.91 

9.75 

8.03 

5. 46 

3.32 


TABLE III. - COEFFICIENT FOR MAXIMUM POSSIBLE ATTENUATION, 
SPINNING MODES, AND CIRCULAR DUCT 



Lobe pattern number, m 

0 

1 

3 

5 

7 

10 

20 

50 

Attenuation coefficient, D 

m 

21. 51 

35. 88 

66. 83 

97.87 

134.9 

184.9 

372.3 

1099 
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